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ABSTRACT 

Relativistic charged particles (CR for cosmic-rays) produced by Supernova explosion of the first 
generation of massive stars that are responsible for the re-ionization of the universe escape into the 
intergalactic medium, carrying an electric current. Charge imbalance and induction give rise to a 
return current, jt, carried by the cold thermal plasma which tends to cancel the CR current. The 
electric field, E = rjjt, required to draw the collisional return current opposes the outflow of low energy 
cosmic rays and ohmically heats the cold plasma. Owing to inhomogeneities in the resistivity, ri{T), 
caused by structure in the temperature, T, of the intergalactic plasma, the electric field possesses a 
rotational component which sustains Faraday's induction. It is found that magnetic field is robustly 
generated throughout intergalactic space at rate of 10~^^ — 10^^^ Gauss/Gyr, until the temperature of 
the intergalactic medium is raised by cosmic reionization. The magnetic field may seed the subsequent 
growth of magnetic fields in the intergalactic environment. 

Subject headings: plasma — large-scale structure of universe — magnetic fields:theory 



1. INTRODUCTION 

Massive stars characterized by high emission of ion- 
izing UV photons are the main co ntributors to the 
proce ss of cosmic re-ionization (e.g., iCiardi fc Ferraral 
|2005[) . At the end of their life they explode as 
Supernovae acc elerati n g, large amounts o f CR pro- 
tons (Krvmsk-^ HQTTI: lAxford et al.l [TOTTI iBeil ITOTI 
IBlandford fc OstrikeiiUQTSH . Since dynamo generat ion of 
magn etic fields operates very quickly within stars (|Reesl 
l2006f) . stars can magnetize their surroundings through 
their own magnetized stellar winds. Diffusive shock ac- 
celeration around Supernovae can, therefore, take place 
even without a pre-existing galactic magnetic field. Ow- 
ing to their much higher energy and diffusive mean 
free path compared to thermal particles, CRs eventu- 
ally escape from the parent galaxy into the intergalactic 
medium. Their escape may or may not be collimated 
by the galaxy disk (e.g. by the magnetic field in the 
wind around stars or a pre-existing galactic field if it ex- 
ists). In any case, the CR protons carry a small but 
important electric current jc which couples them to the 
thermal intergalactic medium and induce magnetic fields 
generation there. The phenomena described in this pa- 
per depend on the electrical current carried by the CR. 
Most of the current is carried at GeV energies, so our dis- 
cussion does not depend upon CR acceleration to TeV or 
PeV energies. Also we neglect CR electrons as they are 
typically two o rders of magnitude less numerous than 
protons (jSchlick ciscr 2002). 

The Larmor radius of a proton of momentum p in a 
magnetic field B is rg = {p/mc){B/10~^^G)~^kpc. In 
our discussion we will consider fields up to lO^^^G and 
protons energies up to IGeV, so magnetic field will place 
very little restriction on intergalactic CR propagation. 
Even a spatial diffusivity unrealistically as small as the 



Bohm diffusivity, in which the mean free path is set equal 
to the Larmor radius, would allow GeV protons to diffuse 
a several Mpc in a 10^ Gauss field after 1 Gyr. Inter- 
particle collisions also place no restriction on propagation 
since the Coulomb collision time for GeV protons in a 
characteristic plasma density ~ 10~'*cm~'^ is ~ lO'^Gyr. 
CR propagation within galaxies is not necessarily as clear 
cut since it is possible that the magnetic field, although 
largely unknown, may conceivably be much larger than 
in the intergalactic medium as result of processes associ- 
ated with stars and supernovae. However, CR appear to 
diffuse freely around our own galaxy with present mag- 
netic fields in the /iC range. Even at energies in the GeV 
range, the presumption is that the galactic CR spectrum 
is determined by a balance between CR production by 
SNR and escape from the galaxy (jHillas 2005 ). More- 
over, there is relatively direct evidence that even in the 
large magnetic field close to the center of our galaxy CR 
diffusion, far exceeding Boh m diffusion, is relatively free 
(|Dimitrakoudis et al.l 120091 ). Hence we base our discus- 
sion on the assumption that CR produced by early su- 
pernovae escape freely from their host galaxy. Once GeV 
protons escape into the intergalactic medium their prop- 
agation is restricted by neither collisions nor magnetic 
field. 

In this paper we show that the electric currents car- 
ried by CRs protons produced in the first generation of 
galaxies that are also responsible for re-ionization of the 
universe is sufficiently large to generate significant mag- 
netic fields throughout intergalactic space. The magnetic 
field is generated by the curl of the electric field associ- 
ated with the return currents induced by the CR parti- 
cle propagation. This process operates most efficiently 
while the intergalactic medium is cold, and effectively 
shuts down once reionization raises the intergalactic gas 
temperature above 10^ K, i.e. when the universe was 
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about 1 Gyr old ()Ciardi fc Ferraral |2005[ ). The gener- 
ated magnetic field is typically of order lO^^^G at the 
end of reionization. This magnetic field may ev entually 
be further enhanced by a turbulent dynamo (Rv u et aLl 
|2008), although CR currents may play a further ro le in 
magnetic fi eld amplification through Lorentz's force (jBelll 

[2onl[2nn5h . 

Magnetic fields are observed in most astrophysical bod- 
ies from planetary scales, to stars and galaxies, and up 
to the largest structures in the universe fZel'dovich et alj 
[1983: Kulsrud fc Zwcibc l 2008 ). Nearby galaxies have 
been known to be magnetized for some time, but re- 
cent studies show that they acquire their fi elds when the 
universe was less than half its present age (jBernet et alJ 
[200l IKronberg et all [2008h . Since the growth of large 
scale magnetic field in a dynamo model is exponential in 
time, this poses serious restrictions on the timescale on 
which galactic-dynamo must operate. 

Evidence for intergalactic magnetic fields is provided 
by Faraday rotation measure and diffuse synchrotron ra- 
diation which reveal the existence of fiG str ong magnetic 
fields in the hot pla sma of galaxy clusters (jClarke et al.l 
120011: iCarilh fc Taylo r 2002). These probes are consid- 
erably less sensitive to fields in adjacent structures such 
as small groups and filaments of galaxies. Nevertheless, 
these fields are expected to be there and large efforts 
are being made in an attempt to measure them. In 
addition, intergalactic magnetic fields affect the prop- 
agation of ultra high energy cosmic ray particles, in- 
troducing potentially measurable effects on th eir energy 
spect r um, arrival directi o n and composition (ISigl et al.l 
[200l IDolag et al.l [200l IHoooer fc Taylml [20m fT 
nally, the presence of magnetic fields in voids can af- 
fect the observed spectrum of extragalactic TeV 7-ray 
sources. Multi-TeV photons are absorbed by the dif- 
fuse extragalactic background light and converted into 
pairs which emit secondary cascade multi-GeV 7- 
rays by inverse Compton on the cosmic microwave back- 
ground. The observed flux of pair produced GeV pho- 
tons can be suppressed if the pairs are deflected from 
the line of si ght by a magnetic field within an energy 
loss distance ( Aharonian et al.]|1994l : lNeronov fc Semikozl 
I2009D . Based on these ideas, in very recent studied! 
the 7-ray spectra of blazars were used to set lower 
limits to the value of magnetic fields in voids in the 
range 5 x IQ-^^ _ 5 ^ iQ-i^Q (|N eronov fc Vovk|[20T0l: 
iTavecchio eraL|[20Tol ). 

The presence of magnetic fields in astrophysical plasma 
requires a mechanism for their generation because, given 
the high mobility of the charges, sustained electrostatic 
fields are non-trivial to set up. However, even a weak 
seed may suffice because magnetic field amplification can 
occur at the expense of the plasma motions through the 
induced electric field, E = —v x B. Specific lower limits 
on the required seed strength depend on the astrophysi- 
cal systems. 

There already exist various scenarios for the gen- 
eration of magnetic fields each with its strengths 
and weaknesses, among others Weibel's instability 
at sh ocks (Mcdvc dev et al.l 120061: iSchlickeiser fc Shuklal 
120031 [Medvedevi (20071); battery effects, e.g. due to 

^ these results appeared after the first submission of this paper 



Compton drag (|Harrisonl 119701: llchiki et al.l I2006D or 
Biermann's me chanism either at cosmic shocks or ion- 
ization fro nts (iSubramanian et all 119941 : iKulsrud et al 



1997: Gnod in et al. 20001) : galactic winds (Bertonc et al 
200a JDonnert et alj 120091) , processes in the early uni- 
verse (j Widrow 20021). In fact, the magnetization of cos- 
mic space is complex and various processes are likely to 
contribute to it, though to a different extent in differ- 
ent environments. One of the strengths of the model 
we describe below is that it generates relatively strong 
macroscopic fields throughout cosmic space, while having 
a simple description based on well understood physical 
processes. 

The remainder of this paper is organized as follows: In 
Section [2] we introduce the resistive mechanism for the 
generation of magnetic field. In Section [3] we describe 
a cosmological simulation to determine some of the pa- 
rameters describing the intergalactic medium that are 
required to quantify the generation of magnetic field. In 
Sections |4] and [5] we discuss the generation of magnetic 
field around individual galaxies using a numerical and 
slightly simplified analytical solutions, respectively. The 
solution is extended to the case of the intergalactic space 
in [5] Finally, we briefly summarize the main flndings of 
this paper in [71 

In the following we use SI units and, unless explicitly 
stated, lengths are expressed in physical units, not co- 
moving units. 

2. RESISTIVE MECHANISM 

Though, as discussed in the previous section, CR prop- 
agation in intergalactic space is unaffected by Coulomb 
collision and magnetic flelds deflections, CRs are not 
completely decoupled from intergalactic plasma. This 
is due to the small current, jc, carried by the CR flux, 
which is assumed mainly to consist of protons. In fact, 
to maintain quasi-neutrality the CR current must be bal- 
anced by a return current carried by the thermal plasma 
jt- While the CR are collisionless with very long mean 
free paths, the thermal particles, due to their low en- 
ergy, have mean free paths shorter than scale lengths of 
interest here. For example, at redshift z ~ 10, when 
the intergalactic plasma had a temperature ~ IK and 
density ~ 10~''cm~'^, thermal electrons had a Coulomb 
mean free path of only 10^V~^ii^ where ij^ is the frac- 
tional ionization. Consequently an electric field in the 
rest frame of the plasma (indicated by a /) is required to 
draw the return current, namely 



E' = VJt 



(1) 



where rj is the plasma resistivity. Since Coulomb col- 
lisions dominate over charge-neutral collisions the re- 
sistivity takes approximately the Spitzer value, 77 = 
65T~'^/^logA flm, where T is the temperature of the am- 
bient thermal plasma in K, and logA ~ 20 is the Coulomb 
logarithm. The Spitzer resistivity is independent of elec- 
tron number density and consequently independent of 
the degree of ionization. 

The electric field transfers energy from the CR flux to 
the thermal plasma through ohmic heating, described by 
the equation 



3 , dT 
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where ks is the Boltzmann constant. Since the charge- 
neutron cohision time is much smaller than the expansion 
time of the universe, in Eq. ([2]) n is the total number 
particles in the plasma, to include heating of neutral as 
well as ionized particles. As a result of ohmic heating 
the temperature of the initially cold thermal plasma is 
raised. 

The electric field also opposes the current carried by 
the CRs. In particular, CRs escape the parent galaxy to 
a distance R only if their kinetic energy at source exceeds 

the potential ai R, (j) = \E\dR. This can determine 
the minimum energy of CR reaching R, giving 



S(^(i?) 



1 



1/2 



(3) 



As will be seen in figure 3, the electric potential is not 
sufficient to limit the escape of mildly relativistic protons 
from the galaxy for the cases considered here, largely 
because we assume that gas within lOOkpc of a galaxy is 
hot (see below) , and therefore has a low resistivity which 
results in a low electric field even though the CR current 
is relatively large close to the galaxy. 

In a uniform medium, the electric field drawing the re- 
turn current jt cannot have a curl and therefore there is 
no generation of magnetic field. If the medium is non- 
uniform on a given scale, the forward and return currents 
can become separated on that scale producing a current 
loop which in turn supports a magnetic field. However, 
the separation of forward and return currents, and the 
corresponding magnetic field, is opposed by inductive ef- 
fects, and the actual growth of magnetic field must be de- 
termined by self-consistent solution of the Maxwell equa- 
tions 



dB 



dt 



= -V X E. 



(4) 
(5) 



The displacement current can be neglected because 
timescales are much longer than the light transit time 
and the system evolves through a series of quasi-neutral 
steady states. Using Ampere's law (|4]) to eliminate the 
thermal current we find that the electric field ([1]) in the 
frame in which the the plasma moves at velocity v, is 

E = —vx B + {ri/ fj,Q)V X B -rjjc- The curl of this electric 
field then produces growth of magnetic field according to 
Faraday's la-«EI 

^ = ^ X {vx B)-V X ( — V X B ) + V X {r]jc). (6) 
ot V Mo / 

The first term on the right hand side transports the 
frozen-in field with the plasma and can stretch and am- 
plify an already existing magnetic field. The second term 
represents resistive diffusion and can easily be verified 
to be insignificant for kpc distances and Gyr timescales. 
Crucially, in contrast, the final resistive term of the equa- 
tion produces magnetic field in a previously unmagne- 
tized plasma. The magnetic field grows wherever the 
resistivity varies perpendicularly to the CR current. A 

^ Adiabatic losses due to cosmic expansion are negligible during 
the ~ 1 Gyr in which the resistive process is efficient. 



larger electric field is needed to draw the return current 
where the resistivity is higher, so the electric field has 
a curl and magnetic field grows. These equations are 
well known in laser-produced plasmas where the current 
is carried by energetic laser-produced electrons in place 
of Su pernovae-produced CRs protons ((Bell fc Kinghaml 
I2OOI . 

Temperature inhomogeneities, resulting in variations 
in conductivity, are naturally present throughout inter- 
galactic plasma due to the growth of cosmological struc- 
ture. In addition, because the volume heat capacity 
is proportional to density, where ohmic heating is sig- 
nificant, inhomogeneous temperature enhancements will 
arise from density inhomogeneities. On the other hand, 
the CR current is expected to be approximately uniform 
because as discussed previously sections, CR are rela- 
tively undeflected by collisions or magnetic fields in the 
intergalactic medium. Thus the magnetic field may be 
generated by the resistive term at a rate 



B 



(7) 



where, Lt = T/\\7T\, is the characteristic temperature 
scale. 

It is clear from the above expression that the resis- 
tive process depends sensitively on the plasma tempera- 
ture through 77. In particular, it operates efficiently while 
the intergalactic medium is cold, and it effectively shuts 
down once reionization raises the intergalactic gas tem- 
perature above 10 ^ K. For a uni verse that is reionized 
at redshift z > 6 (jCiardi fc Ferrara 2005.). the available 
time is of order a Gyr. Note that close to the star-forming 
galaxies the gas temperature is > lO^'K due to both pho- 
toheating from the ionizing flux escaping the galaxies and 
the ohmic heating. For the above choice of the reioniza- 
tion epoch, the typical size of HII regions at z 10 is 
about 100 kpc (or 1 comoving Mpc. iZahn et al.ll2007h . 

In order to quantify the importance of the resistive pro- 
cess in the following we first use a numerical simulation 
of structure formation to compute the scale of the tem- 
perature variations in the intergalactic medium. We then 
estimate the generation of magnetic field around individ- 
ual star-forming galaxies at the epoch of re- ionization, in- 
cluding self-consistently the effect of the induced electric 
fields on the escaping CR particles and ohmic heating. 
Finally, we estimate the distribution of magnetic fields 
produced in intergalactic space. 

3. COSMOLOGICAL SIMULATION 

We extract characteristic scale-length Lt and the char- 
acteristic range of values taken by Lt from a cosmologi- 
cal simulation of structure formation which includes hy- 
drodynamics, the relevant thermodynamic processes for 
the diffuse baryonic gas, dark matter and gravity. The 
simulation does not include the cosmic ray effects dis- 
cussed in this paper. The simulation uses a directionally 
un-split higher order Godunov's method for the hydro- 
dynamics, a time centered modified symplectic scheme 
for the collisionless dark matter and we solve Pois- 
son's equation with a second o rder accurate discretiza- 
tion (jMiniati fc Colellal l2007aD. Atomic cooling, with 
rates as given in IHui fc Gnedh] ([l99%_ is included ac- 
cording to the algorithm proposed in ( Miniati fc Colellal 
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Fig. 1. — Baryonic gas density distribution. 

l20071ilMiniatil[20l0l) . although prior to reionization adi- 
abatic expansion of the universe is the most relevant 
cooling process. We adopt a fiat ACDM universe with 
the following parameters normalized to the critical value 
for closure: total mass density, ftm — 0.2792, bary- 
onic mass density, flf, — .0462, vacuum energy density, 
= I - n.m = 0.7208 ()Komatsu et al.|[2009t ). In ad- 
dition, the normalized Hubble constant \s h = TJo/lOO 
km s^^ Mpc^^ = 0.701, the spectral index of primordial 
perturbation, ria — 0.96, and the rms linear density fluc- 
tuation within a sphere with a cqmovi ng radius of 8 h^^ 
Mpc, (78 = 0.817 (|Komatsu et al.ll2009( ). W e generate the 
initia l conditions with the grafic2 package (|Bertschingei] 
l2001h . We use a computational box of comoving size 
L = \h~^ Mpc, discretized with 512'^ comoving compu- 
tational cells, providing a nominal spatial resolution of 
2/i~^ comoving kpc for the field components, and 512^ 
particles with mass 4.8 x \Q^h~^ Mq for the collisionless 
dark matter component. At redshift z = 6 the box size in 
physical units is ~ 143 /i~^kpc and, likewise, the nominal 
spatial resolution is ~ 285 /i~^pc. Thus, structures with 
spatial scales between 1-10 kpc should be adequately re- 
solved for our purposes. 

Figure [1] shows the distribution of the baryonic gas 
density at z ~ 10 on a two dimensional slice across the 
simulation box. One can recognize a few high density 
collapsed structures where gas is rapidly cooling. These 
are the sites where stars and galaxies form and CR are 
eventually produced. However, most of the gas (~ 97% 
by mass) is still in the diffuse phase, with a density within 
a factor a few of the mean value. Figure [5] shows the 
occurrence of spatial scale-lengths in the range 1-10 kpc 
in the low density, low temperature gas which occupies 




Fig. 2. — Temperature (top) and its characteristic length scale 
(bottom) as a function gas density in units of its average value. The 
vertical bars represent the asymmetric root-mean-squared fluctua- 
tions about the mean. 

most of the intergalactic space. 

Note that the residual fraction of free electrons in 
the intergalactic medium after recombination keeps the 
gas temperature T locked to the CMB temperature, 
TcMB = 2.725(1 -I- z) K, until red shift 1 -{- z = 
142(f76/iV0.024)2/5 ~ 140 (|Peebleslll99l . After that the 
temperature of the diffuse intergalactic medium drops 
due to adiabatic expansion following the 7— law for a 
monoatomic gas with 7 — 5/3. 

4. SOLUTION AROUND INDIVIDUAL GALAXIES 

The growth of magnetic field around an individual 
galaxy is described by Eq. ([7|). To solve this equation 
in the following we use the temperature scale-length de- 
rived from the simulation in the previous section as a 
reference value for Lt- In addition, we will determine 
the CR current produced by a galaxy of a given lumi- 
nosity. The current determines the electric field which 
sets the minimum momentum of the escaping CRs, Eq. 
([3]), from which the current itself depends. In addition, 
the electric field is effectively related nonlinearly to the 
current because the resistivity is affected by ohmic heat- 
ing, Eq. ([2]). So all these equations need to be solved 
self-consistently. 

The luminosity of a typical bright galaxy at redshift 
z > 6, is ~ 5.2 X lO^^WHz"^, corresponding to a 
star formation rate Af* = 6.5AfQyr~-'^ (|Bouwens et al.l 
|2007| ). If we conservatively assume a Salpeter ini- 
tial mass function, the energy released by SN explo- 
sions per unit mass of formed stars is, Esn — 5.4 x 
10*2 J Afg^ A fraction, ~ 30%, of this energy is 
typically converted into CR particles, implying a CR- 
luminosity, ~ 3.2 x lO^^ W (ec/0.3) (L/L*), for a 
galaxy of luminosity L. The corresponding CR en- 
ergy flux at a distance R from the galaxy is Qc — 
2.7 X 10-5(ec/0.3)(L/L*)(i?/kpc)"^Wm-^ where we 
have idealized the galaxy as an isotropic point source 
of CR. If the CRs are coUimated within a solid angle fic 
then the CR current intensity is higher by a factor Air/Q.c, 
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although this does not change the quahtative picture. 
The CR momentum distribution can be expected to be 
p~^-^ as typically observed for CR production by shocks 
in the Galaxy. With this power-law, most of the energy 
resides in mildly relativistic protons. The electric current 
is predominantly carried by the lowest momentum CR. 
In the following analysis, the electric field is allowed to 
inhibit CR propagation, although we find that inhibition 
is negligible for our conditions. If only CR with mo- 
mentum greater than Pmm reach a radius R, the current 
carried by CR at that radius is approximately given by 
jc = 5 X lO~^°QciPmtn/mpc)~° y{pmin/mpC+ 1) Amp 
m~^ where irip is the proton mass. Taking the efficiency 
of CR production to be fixed at Cc = 0.3 



J, ~ 5.3 X 10-^^ ( ^ 



1 + 



R 

kpc 

Prnin 



Pn 



niriC 



-0.3 



Amp m (8) 



This CR current must be balanced by an equal but op- 
posite return current drawn by an electric field E = rjjc, 
giving 



~ 6.9 X 10"" 



i- 

Vkpc 



Pn 











Vtnin \ 




nipC J 





-3/2 



-0.3 



V m" 



(9) 



So we can now solve numerically the coupled system of 
equations ([2|), ([3|), dH), ([9]), to determine the propagation 
of CR through the intergalactic medium surrounding a 
star forming galaxy, the associated electric field and the 
ensuing generation of magnetic field. Our 'standard' cal- 
culation is as follows. We treat the galaxy as a spherically 
symmetric source of CR which is constant in time. We 
set an inner radial boundary to the computational grid 
at rmin = lOkpc which corresponds to the notional ra- 
dius of the galaxy. We assume that all CR produced by 
the galaxy escape through this boundary for the reasons 
given in Section I. As will become clear below, changing 
Tmin has vcry little effect on the results for the inter- 
galactic medium on Mpc scales. The CR flux is specified 
at the inner boundary as a power law (oc p~^'^) with 
a minimum momentum of O.lmpC which corresponds to 
a proton velocity a few times the shock velocity of an 
expanding young supernova remnant. We set the initial 
intergalactic temperature to 2 x lO'^K within a radius 
of lOOkpc to represent heating due to ionization by the 
galaxy. Outside this radius we set the initial temperature 
to 1 K, with the temperature changing between the two 
regions over a distance of lOkpc. We set the intergalactic 
density to 10~^cm~^. 

4.1. CR propagation and magnetic field generation 

Figure 3a shows the results for our standard calcula- 
tion, as defined above, for CR propagation and plasma 
heating as a function of radius for a galaxy with typical 
luminosity L = at time t = 1 Gyr. The horizon- 
tal scale is logarithmic and stretches from the edge of 
the galaxy at lOkpc to a distance of lOMpc which is 
characteristic of the distance between bright galaxies at 
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Fig. 3. — Spatial profiles at radii between lOkpc and lOMpc. a) 
standard calculation, b)-f ) as for the standard calculation with one 
parameter changed as indicated. The units are CR current den- 
sity j in 10~^*Am~^, temperature T in degrees K, magnetic field 
B in 10-22T (10-l*G), electric field E in lO"!* Vm"!, electric 
potential </> in Volt. 



this epoch. Ohmic heating is effective out to a radius 
of ^ 3Mpc where the electric and magnetic fields reach 
their maximum values. The magnetic field is less within 
this radius because ohmic heating reduces the resistiv- 
ity and the growth of magnetic field. Within lOOkpc 
the magnetic field is much smaller because the temper- 
ature is initialized to a much larger value to represent 
ionization by the galaxy. The magnetic field is largest 
at 3Mpc, but the maximum is broad, and fields above 
10" ^'^G extend from lOOkpc to lOMpc. Figure 3b dis- 
plays the profiles at the earlier time of O.lGyr for the 
same parameters. The fields are slightly reduced but 
still reach 2.5 x 10~^''G, once again with a broad max- 
imum. This shows that most of the field generation oc- 
curs quickly before ohmic heating reduces the resistivity. 
Fields in the range 10~^^ — 10~^^G are produced even if 
suitable conditions for growth apply for much less than 
IGyr. 

We take figure 3a as the standard calculation and then 
vary individual parameters in figures 3c-3f. Fields of a 
similar magnitude are produced if the initial tempera- 
ture is raised to lOK (figure 3c), the CR luminosity of 
the galaxy is reduced tenfold (figure 3d), the density is 
increased to 10^'^cm"'^ (figure 3e), or the minimum CR 
energy is increased to Pmin = nipC (figure 3f). Clearly the 
production of fields in the range 10"^^— 10~^^G is robust. 
The results indicate that fields as large as ^ 10~^°Tesla 
(10~^®Gauss) are produced in a few to several ^Mpc'^ 
volume surrounding galaxies if the temperature or den- 
sity vary on kpc scales. 

The plots of the electric potential (j> in figure 3 show 
that the electric field does not inhibit the escape of CR 
from the galaxy. 

5. ANALYTIC SOLUTION FOR A FIXED CURRENT 

For the cases considered above, the electric potential 
does not inhibit CR transport. If the galaxy can be as- 
sumed to produce CR at a constant rate, the CR current 
is constant at any point outside the galaxy. The evolu- 
tion of the temperature and magnetic field at that point 
is then determined by the equations 
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dt Lt ' dt 3nki 
The solution is 




(10) 



(11) 



where Ti is the initial temperature in Kelvin, Bi = 



{r]inkBTit/L^) 



2 ^1/2 



Ji 



{nkBTi/rjity/^ , and ryi is the 



resistivity at temperature Ti. The magnetic field is 
largest when jc = 5.15ji. When jc exceeds 5.15ji, ohmic 
heating increases the temperature, and reduces the resis- 
tivity, which in turn reduces the growth of magnetic field. 
Ohmic heating is marginally important when jc — 5.15ji. 
At this CR current the magnetic field achieves a maxi- 
mum value, Bmax — 1.046-Bi. Expressed in more mean- 
ingful terms, after a time t in Gyr, the maximum field 
is 



Anax = 8.2 X 10-" 



10-4cm-3 



1/2 / 



VGyr 



K 

1/2 



-1/4 



(12) 



G. 



The maximum field occurs where the CR current (jc 
5.15ji) is 



Jmax — 3 X 10 



-20 



10-W^'' ' ~ ' ^^^^ 



K 



-1/2 



Ar 



and the maximum field occurs at a distance -Rmax from 
a galaxy with luminosity L, where 



R„ 



1.9 



10-4cm-3 



-1/4 /Ti\-^/^ / L 



0.1m„c 



-0.15 



O.lTOpC 



K 



-1/2 



(14) 



1/4 



(g^) "^'"^ 



in a medium with a proton density n. With the proviso 
that i?max lies in a realistic range, Bmax is independent 
of the emissivity of the CR source, the CR momentum 
Pmin and the fraction of CR escaping the galaxy. i?max 
is proportional to the inverse of the temperature scale- 
length, but all other dependences are relatively weak. 
The weak dependence on the initial temperature implies 
that significant magnetic field is generated even if the 
initial temperature is greater than IK. Hence, the char- 
acteristic magnetic field is robustly of the order of 10~^^ 
to 10~^^G when the effect of Ohmic heating is included 
in the calculation. The maximum magnetic fields and 
their spatial locations in figure 3 obey these equations. 
The weak dependencies on various parameters results in 
the broad extent of the maxima in figure 3. 





Fig. 4. — a) Distribution functions of current and magnetic 
field for galaxies randomly distributed in space according to the 
Scliechter function. Apart from the luminosity the Monte Carlo 
calculation adopts the standard parameters used in figure 3a. b) 
distribution functions when the CR luminosity of all galaxies is 
reduced tenfold. 

6. SOLUTION FOR THE INTERGALACTIC MEDIUM 

The number density of star forming galaxies at redshift 
z > 6, with luminosity, L, and per luminosity interval, 
dL/L^,, is well described by a Schechter function 



$ (i) - ( ^ 



(15) 



with the following parameters: ~ 5.2 x 

10^^ergs-^Hz~\ - IQ-^Mpc" "^ a ~ 1.77 

(iBouwens et al.l 120071: lOesch et al.l 120091: iBouwens et"all 
120101 ). is the typical luminosity of a bright galaxy 
and roughly corresponds to their number density, a 
is the slope of the distribution at the faint luminosity 
end. 

The electric current at any given point in intergalac- 
tic space is contributed by galaxies within the cur- 
rent horizon, Rjh- As long as the diffusion coeffi- 
cient remains much larger than Bohm's value, as ob- 
served in the much more turbule nt interstellar medium 
of the Galaxy dProge fc Kartavv kh 2009), even for B ^ 
10~^^G, Rjh is of order a few tens of Mpc, i.e. larger than 
the average distance between bright galaxies. Magnetic 
field generation is thus expected to be dominated by the 
nearest luminous galaxy. This is the case even if the cur- 
rents are beamed with an opening angle 9 ~ 0.5 rad. In 
this case, the number of CR current sources visible from 
any given point is. 



Nr. 



^,T{2 ~ a)R%{9y47T) ^ 1, 



(16) 



for Rjh > a few x 10 Mpc. Therefore, each point 
in space is exposed to the CR current from about one 
galaxy. Note that because the faint-end slope of the 
above Schechter's function is steep, i.e. a < 2, the CR 
output per luminosity log-interval, ecL,(L/L*)2$(L) ~ 
eci*(.Zv/L*)2~"$,, is approximately constant. However, 
the average distance between faint galaxies of luminosity 

(dL) = [i<f (i)]"'/' oc l("-i)/3, (17) 
so that their CR current at this distance scales as 

jc(L,dL) cxe,i/(dL)2 aLi/2^ (18) 

i.e. fainter galaxies are slightly less efficient at magnetiz- 
ing their surroundings. 

We now make a more quantitative estimate through a 
Monte-Carlo simulation in which galaxies are distributed 
randomly in space according to the Schechter function as 
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given in equation (jlSp . The currents from each galaxy 
are added vectorily and the calculation repeated many 
times to build-up a distribution function for the CR cur- 
rents found at any random point in space. The magnetic 
field is then calculated using equation 11. Figure 4a gives 
the distribution functions for the current and the mag- 
netic field for the standard parameters assumed in figure 
3a. The distribution functions are defined as the number 
per logarithmic interval in j or B. The magnetic field 
calculation assumes that the initial temperature is IK 
everywhere. It neglects heating during ionization, which 
was included in figure 3 within a radius of lOOkpc, but 
this neglect has little effect on the distribution function, 
firstly because the ionized volume is relatively small com- 
pared with the characteristic distance of lOMpc between 
luminous galaxies, and secondly because, even starting 
from a low temperature, ohmic heating reduces the re- 
sistivity, and therefore the magnetic field, in regions close 
to galaxies. For fixed parameters such as the tempera- 
ture scale-length Lt and the plasma density n, the mag- 
netic field cannot exceed a maximum value as given in 
equation 12. This produces an abrupt cut-off to the dis- 
tribution function. The spike in the distribution function 
at the cut-off in figure 4a reflects the broad spatial extent 
of the maxima in figure 3a. The magnetic field is close 
to the maximum value over a wide range of distances 
from individual galaxies. In reality, as shown in figure 2, 
the temperature scale-length Lt takes a range of local 
values, and this will smear the cut-off. The effect of de- 
creasing (increasing) the temperature scale-length Lt by 
a factor A would be to move the distribution function for 
B a corresponding factor A to the right (left) in figure 4a, 
while leaving the distribution function for j unchanged. 

In figure 3 it can be seen that the value of the maximum 
magnetic field is insensitive to changes in various param- 
eters. However the position of the maximum moves in 
radius. For example, the maximum magnetic field oc- 
curs at IMpc when the CR luminosity of all galaxies is 
reduced by a factor ten (see figure 3d where L = O.IL*). 
The effect of this on the distribution function for B can 



be seen in figure 4b. When the luminosity is reduced, 
there are still regions of space in which the magnetic 
field reaches the maximum value of 8 X 10" ^'^G, but the 
field is characteristically 10~^^G throughout most of the 
volume between galaxies. Similarly, the effect on the dis- 
tribution function of changing other parameters can be 
deduced from figure 3 in conjunction with equation 12. 
Equation 12 gives a good estimate of the maximum mag- 
netic field, while figure 3 indicates the volume of space 
that is filled by fields close to the maximum value. 

7. SUMMARY 

In conclusion, we have shown that cosmic rays pro- 
duced by supernovae at the epoch of reionization when 
the first stars and galaxies form can be expected to have 
a substantial impact on the intergalactic medium out to 
distances of a several Mpc from the galaxy. Through the 
return electric current carried by the thermal plasma, 
the medium is heated ohmically and magnetic field is 
generated in the range 10~" to lO-i^Gauss. Because of 
self-compensation through ohmic heating, the maximum 
magnetic fields are relatively insensitive to parameters 
such as the galaxy luminosity, the efficiency of CR pro- 
duction, and the ambient intergalactic temperature and 
density. The magnitude of the magnetic field is however 
sensitive to the scale-length Lt of temperature variations 
in the intergalactic medium. Cosmological simulations 
support an assumption of scale-lengths in the range 1- 
lOkpc. Self-consistent cosmological simulations including 
cosmic ray effects are a necessary next step in develop- 
ing this theory, but our argument extends the possibility 
that resistive magnetic field generation, driven by cosmic 
ray streaming, may generate seed magnetic fields which 
might subsequently be turbulently amplified to form the 
intergalactic magnetic field at the present epoch. 



FM acknowledges several useful discussions with Dr. 
P. Oesch. 
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